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ABSTRACT 

The Sunyaev-Zel'dovich (SZ) effect has long been identified as one of the most important 
secondary effects of the Cosmic Microwave Background (CMB). On the one hand, it is a 
potentially very powerful cosmological probe providing us with additional constraints and 
on the other hand it represents the major source of secondary fluctuations at small angular 
scales (£ > 1000). We investigate the effects of the SZ modelling in the determination of 
the cosmological parameters. We explore the consequences of the SZ power spectrum com- 
putation by comparing three increasingly complex modelling, from a fixed template with an 
amplitude factor to a calculation including the full cosmological parameter dependency. We 
also examine the dependency of the cosmological parameter estimation on the intra-cluster 
gas description used to calculate the SZ spectrum. We show that methods assuming an SZ 
template bias the cosmological parameters (by up to 2<x on <x 8 ) when the cosmology used in 
the template deviates from the reference one. A joint CMB-SZ analysis with a full cosmolog- 
ical dependency of the SZ spectrum does not suffer from such biases and moreover improves 
the confidence intervals of cr% and Q.& m h 2 (2.5 and 2 times respectively) with respect to a pure 
CMB analysis. However, the latter method is quite sensitive to the intra-cluster gas parameters 
and hence requires extra information on the clusters to alleviate the induced biases. 

Key words: cosmology: theory - methods: statistical - cosmic microwave background - 
cosmological parameters - galaxies: clusters: general 



1 INTRODUCTION 

Over the last decade, the low multipole observations of the Cos- 
mic Microwave Background (CMB) angular power spectrum, in 
intensity and polarisation, show that a ACDM concordance model 

J II "T 

describes accurately the Universe (Komatsu et al. 2009) with the 
basic cosmological parameters constrained dDunklev et al.ll2009h 
with a precision of the order of a percent. To obtain even bet- 
ter constraints on "standard" cosmological parameters and in or- 
der to further constrain the cosmological model (dark energy, 
running of the spectral index, etc.), various experiments are col- 
lecting t _data_esrjecmllyayu BIMA, ACBAR , 



SZA) dHolzapfel et alj200cllDawson et alj200ll : |padin et al.l20oTl ; 
iRunvan et alj2003l ; lMuchovei et al.ll2007l) . These experiments dis- 
cussed the existence of a power excess at large t values that could 
be accounted for by point sources, SZ effect or more exotic physics 
-non standard inflation, primordial voids, feat ures in the primordial 



spectrum, primordial non- Gaussianity, etc- dElgar0v et al. 



2002; 



Cooray & Melchiorri 2002; Griffiths et al. 2003; Bond et al. 2005; 
Dawson et alj 120061; IDouspis et al.1 l2006l : iReichardt et alj 120091 ; 
Sievers et al . 2009; S harp et alj|2009h . This highlights the neces- 



sity to carry out a consistent analysis of the CMB signal, based on 
a model that describes both primary anisotropics and also the sec- 
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ondaries arising from the interaction of CMB photons with matter 
between the last scattering surface and the observer. 

In the case of experiments that have access to high multi- 
poles (C > 1000), the secondary SZ anisotropies and the point 
sources contribution will dominate over the primary CMB. A 
joint analysis of the CMB and SZ power spectra is thus nec- 
essary. Different parameterisations of the SZ spectrum have al- 
ready been used to analyse the data from WMAP, CBI, ACBAR, 
and BIMA dSpergel etai] | 2003l; iBond et all 120051; IDouspis et alj 
20061: ISpergel et alJ|2007h Ikuo et alj 120071 ; IReichardt et al.ll2009l 
Dunklev et alj|2009l ; Isievers et alJl2009TT In the present analysis, 
we address the issue of the calculation of the SZ spectrum used 
to fit the data. We examine the effect each SZ description (power 
spectrum and intra-cluster gas model) induces on cosmological pa- 
rameter estimation in terms of accuracy and possible biases. 

In section 2, we briefly introduce the thermal SZ effect and 
the calculation of its power spectrum. We then present in sec- 
tion 3 the different methods used to jointly fit the CMB + SZ 
data. We compare and discuss in section 4 the accuracy of the 
parameter estimation for each of these methods and the possible 
induced biases. We summarize our results in section 5. Through- 
out the study we assume a flat ACDM co smological model an d 
use the WMAP5 cosmological parameters dKomatsu et al.l l2009) : 
o- 8 = 0.817, n s = 0.96, Sl m = 0.279, Q. h = 0.046, h = 0.701. 
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2 THE THERMAL SZ ANGULAR POWER SPECTRUM 



The Sunyaev-Zel'dovich (SZ) dSunvaev & Zel'dovichll 19721) is the 
main secondary anisotropy source at arcminute scales. It is made 
of two terms : the first is the thermal SZ (TSZ) due to the in- 
verse Compton scattering of the CMB photons off the hot electrons 
in the intracluster gas, and the second is the kinetic SZ (KSZ), a 
Doppler shift due to the proper motion of clusters with respect to 
the CMB. The KSZ power spectrum is approximately 2 orders of 
magnitude smaller than the TSZ. We thus neglect its contribution 
in the present analysis. The TSZ has a characteristic spectral sig- 
nature, namely a decrease of the CMB intensity in the Rayleigh 
Jeans part of the spectrum and an increase in the Wien part that 
is due to energy transfer from the hot intracluster electrons to the 
CMB photons. This characteristic frequency signature is given, in 
the non-relativistic approximation, by f(x) = [-t^ry _ 4j, where 

v — — 

In the context of the analysis of CMB data, it is necessary 
to take into account the contribution of the TSZ effect. The most 
useful tool to account for this contribution as a function of the an- 
gular scale on the sky is the TSZ power spectrum. It can be ob- 
tained either from hydrodynamical simulations or from analytical 
calculations. Each of these approaches has its own limitations and 
drawbacks. In the first case, the amplitude of the SZ spectrum as 
well as its shape are sensitive to the simulation characteristics. The 
box-size and resolution, or smoothing length, of the simulation af- 
fect the relative amplitudes on large and small scales respectively. 
When the box-size is too small, the number of massive clusters 
is underestimated and so is the SZ power at large scales. The 
resolution of the simulation, or the smoothing length, artificially 
decreases the SZ power at small angular scales (see for example 
IWhite et al.l 12002). The physical model used to describe the gas 
is an other important source of alteration of the SZ spectrum. For 
example, it has been showed that preheating (due to energy feed- 
back from supernovae for instanc e) as well a s radiative cooling, 
that depends on the gas metalicity dDolag et alj|2005h , respectively 
increase s or decreases the am plitude of the TSZ spectrum by a fac- 
tor of 2 ( Ida Silva et alj2 001). The influence of radiative cooling on 
th e SZ spectrum has al so been studied with an analytical treatment 
by |Zhang&Wuld2003l) and has been shown to significantly affect 
the amplitude of the SZ spectrum. 

In the second case, the analytical calculation of the TSZ 
power spectrum is based on two major ingredients: the halo mass 
function, along with a model for the intra-cluster gas distribu- 
tion within these halos. Neglecting the correlation between haloes 
which is much smaller than the Poisson term or the CMB signal 
dKomatsu~& Kitayama il999h. the TSZ angul ar power spectrum can 
be calculated as in Ko matsu & Seliak U2002h : 



C 



dz— 
o dzdQ. J M 



dM d -^\ M M,zt. (1) 



The mass function n(M,z) is given by the theoretical expression 
of Press & Schechter (1972) or by fitting formulae to N-body 
numerical simulations (Sheth&Torme nll999l ;l Jenkins et al.ll 20011 ; 
IWarren et alj|2006t) . The different mass function can lead to differ- 
ent predictions in terms of the predicted numbe rs of halo, and thus 
may i nduce differences in the power spectrum dKomatsu & Seliakl 
2002). The shape of the TSZ angular power spectrum depends on 
the intra-cluster gas distribution and properties through the two- 
dimensional Fourier transform on the sphere of the 3D radial profile 



of the Compton y-parameter for individual clusters 
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where yso(r) = cr T ^4^n c (r). Several models for the electronic 
distribu tion w e (r) can be used. The most co mmonly used are the fi- 
profile (Cavaliere & Fusco-Femiano 1976) and the polytropic gas 
distribution dKomatsu & Seliakll200lh and extensions of these two. 
The m ean gas mass fraction, / g is taken 0.086/h (LaRo que et al.l 
2006). As for the intra-cluster gas temperature T e , it is usually 
assumed to be isothermal and equal to the virial temperature : 
k\,T e = °^"' pM , where /i is th e mean molecular weight and r vir is cal- 
culated using the formula in Bryan & Normari] dl998l) . The shape of 
the TSZ angular power spectrum also depends on the cosmologi- 
cal model through the comoving volume V c , the angular diameter 
distance D A and the mass function n(M, z). 



y3D(r)- 



Cr/D A 



r~dr, 



(2) 



3 SZ IN CMB ANALYSIS 

The estimation of cosmological parameters from the CMB would 
ideally require a pure primary signal. However the measured CMB 
data will contain additional contributions, of which the TSZ domi- 
nates. The TSZ characteristic frequency signature allows us in prin- 
ciple to remove the secondary contribution of the detected galaxy 
clusters from the CMB in multi-frequency experiments. However, 
on the one hand the residual SZ signal, if not taken into ac- 
count in the analys is, biases the cosmological parameters see (see 
iTaburet et aT]|2009h . and on the other hand taking it into account 
requires a good knowledge of the selection function. An alternative 
method to extracting the SZ contribution and modelling the resid- 
uals is to model the total CMB plus TSZ spectra and to determine 
the cosmological parameters using both primary and secondary sig- 
nals. This is the approach used in all the present high I study. 
We can naturally fit the total signal Cl ot with the expression 

Cf = C™ B 0) + Cf 0) (3) 

with Cf z 0) given by Eq. (TJ as in Douspi s~et alj d2006h , where 6 
stands for the set of cosmological parameters we want to determine. 
The main advantage of this, hereafter method 3, is that it includes 
the full cosmological parameter dependency of the SZ power spec- 
trum. Parameterisations of the SZ power spectrum are also used to 
fit the CMB data. They are very efficient from the point of view of 
computation time since they do not require the full calculation of 
the SZ spectrum at each step of the Monte Carlo Markov Chains 
(MCMC). As a result, parameterisations of the SZ spectrum can 
reduce the chains convergence time by a factor of the order of 2. 
However, the drawback is that they do not reflect the full cosmo- 
logical dependency of the SZ spectrum. Two methods fall in this 
category. In the first one, the total CMB power spectrum can be 
fitted with 



cf = c™ B (d) + A sz c t (e ), 



(4) 



hereafter method 1, where A S z is an amplitude factor multiply- 
ing an SZ spectrum template C[0o) calculated analytically for a 
given cosmology described by the set of cosmological parameters 
8 and intra-cluster gas distribution, or obtained from a given nu- 
merical simulation (i.e. for a given cosmology and gas physics). In 
this analysis the templates we used were those used by the WMAP 
team. 

The last method accounts for the main variations in the SZ 
spectrum amplitude, with the cosmological parameters erg and Sl^h 



3 



following Komatsu & Selia 2 J2002h . In this method 2, the total 
CMB spectrum is fitted with 



CT = C™ B 0) + tr\ (Q b h) 2 C'0 O ), 



(5) 



reference 

WMAP 

- 0=0.8 
T =0.93 x 



where C' ( (6o) is the SZ spectrum for a given cosmology and infra- 
cluster gas distribution. 

In the following (Sect. 14. U . we investigate the effects on the 
cosmological parameter estimation of the three methods introduced 
above. For methods 1 and 2, we used the same SZ spectrum tem- 
plate as that used by the WMAP team Jjp ergel et al] l2007h and 
which is based on iKomatsu & Seliakl bOOll l2002l)~~For method 
3, we computed the SZ spectra using the lSheth & Tormenl fl999r) 
mass function and the spherical yS-profile with p = 2/3 for simplic- 
ity. Furthermore in Sect. 14.21 we study the effects of varying this 
profile on the cosmological parameter estimation. 

In order to compare the three different methods, we create 
mock data, at 100 GHz, in the form of temperature and polarisation 
power spectra containing primary CMB anisotropies and SZ contri- 
bution from all clusters using equation Q] We did not consider the 
polarisation induced by clusters sinc e it is negligible compared to 
the primary one (e.g. lLiu et al.120050 . In this study, we do not anal- 
yse the CMB data produced from multi-frequency observations af- 
ter component separation. This would require to monitor precisely 
the residual signal after component separation which is shown to 
mix all foregrounds and may differ from one component separa- 
tion method to another. Such an approach was followed by J. -A. 
Rubino-Martin and gives similar results as ours (private commu- 
nication). We have rather chosen to use the best channel for CMB 
study : the 100 GHz channel. In this channel, galactic foregrounds 
(free-free, dust emission, synchrotron radiation) will contaminate 
the CMB minimally as well as extragalactic radio and IR point 
sources. We also consider a Planck-like gaussian and uncorrelated 
noise power spectrum with a 9.5 arc minutes beam. We ran M CMC 
analyses, using the CosmoMC code (Lewis & Bridle 2002) with a 
modified version of the CAMB code dLewis et al.ll2000h including 
a module that calculates the SZ power spectrum with its full cosmo- 
logical dependency. We followed for this module the computations 
detailed in Sect. 2. Given the instrumental beam we considered, we 
limited the MCMC analysis to multipoles smaller than ^ max = 3500. 
We used a gaussian likelihood function since the probability distri- 
bution function of the SZ spec trum is well approxima ted by a gaus- 
sian for a Planck-like survey dZhang & Shethll2007|) . The MCMC 
were carried out on the set of parameters 6: Q. b h 2 , Q<imh 2 , the ratio 
of the sound horizon to the angular diameter distance 100 x 6 h , the 
optical depth at reionisation t, the spectral index n s , the CMB nor- 
malisation A 6 , and the SZ normalisation factor Asz when method 1 
is used. The deduced parameters are Q. A , O m , cr 8 , z rc , H and the age 
of the universe. To ensure that the MCMC r uns have converged w e 
used the convergence criterion introduced by Dunklev et al. ( 2005). 



4 RESULTS 

4.1 Effects of the SZ modelling on cosmological parameter 
estimation. 

As a first step we examine the choice of the SZ template (used in 
methods 1 and 2). We thus present the MCMC results, when the 
SZ template deviates from the real SZ angular power spectrum (i.e. 
that used to produce the mock data). To account for this possibility, 
the SZ template we use to fit our simulated data is based on the 
one used by the WMAP team dSpergel et al.! : 2007). We normalized 



Figure 1. Theoretical SZ reference spectrum (fi = 2/3, T e = 7\„>, solid 
black line), SZ WMAP template (dotted blue line) normalised to the t = 
2000 value of the reference spectrum, SZ template of the reference model 
but with a steeper fi distribution of the intra-cluster gas (fi = 0.8, red dot- 
dashed line) and SZ spectrum when the electronic temperature is underesti- 
mated by 7%. The black dotted curves represent the Icr envelopes assuming 
a Planck like experiment at 100 GHz with a 9.5 arcmin beam. 



this template so that its amplitude at £ = 2000 equals the ampli- 
tude of the SZ spectrum employed to create the mock data. At such 
multipoles, the gas properties affect more the SZ amplitude than 
does the cosmology. As a result, this normalisation ensures that the 
main difference between the SZ spectrum used to create the data, 
and the template used to fit them in methods 1 and 2, is due to the 
cosmology. 

We represent in figure [2] the one-dimensional distribution of 
cosmological parameters obtained with the three methods. Using 
the SZ parameterisation of method 1 (long dashed blue lines) en- 
larges the lcr error bars on the parameters and biases the determi- 
nation of some parameters. The values of fld m /i 2 and A s are biased 
by a few tenths in units of the standard deviation. This translates, 
for the derived parameters, into a one-sigma bias on <x 8 , and a few 
tenths in units of the standard deviation for C1 A , fl m and H . These 
biases are due to the fact that method 1, by construction, allows 
to vary only the amplitude of the SZ power spectrum and not its 
shape. 

When employing the SZ parameterisation of method 2, the 
confidence intervals are improved. This can be explained by the 
fact that in method 1 we constrained an extra parameter (A S z) while 
in method 2 the SZ amplitude is given by a combination of <t$ and 
Q. b h. However, as shown in figure[2] the biases on some cosmologi- 
cal parameters become more important. The values of Q. b h 2 , Cla m h 2 
and A s are biased by a few tenths in units of the expected precision 
which, for the derived parameters, leads to a two-sigma bias on <x 8 , 
whereas the values of £1 A , Q. m and H suffer from a one-sigma bias. 
These results can be explained by two facts : first the shape of the 
SZ template in method 2 does not vary, and second its amplitude 
at £ < 2000 is higher than the SZ spectrum used to produce the 
data (Fig. [T). The latter forces the estimated erg to be biased to- 
wards a low value in order to decrease the SZ amplitude (given by 
<r 8 (O b /i) 2 , see Eq. [5]) and thus to fit the data. For the same rea- 
son the Q. b and Hq values are also biased. As for the second point, 
the "fixed" shape of the SZ template, deviating from the reference 
one, forces variations of fi m , S1 A and Sl h so that the primary CMB 
spectrum compensates the shape modifications. 
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Figure 2. One dimensional parameter distribution. The curves represent the distributions for the CMB+SZ signal fitted with the 3 different methods. Long 
dashed blue line : method 1, dot-dashed green line : method 2, solid black line : method 3. The SZ template used in method 1 and 2 (blue dotted line in figure 
[TJ differs from that used to create the data (black line in figure[T}. The vertical red lines represent the input values of the parameters we used to create our 
mock data. Note the substantial biases that can affect the estimation of erg when using method 1 or 2. 



Method 3, which does not require any SZ template with a fixed 
shape, gives unbiased determination of the cosmological parame- 
ters as shown by the black curves in Fig(2] The accuracy obtained 
with method 3 is also better than that obtained with methods 1 and 
2. In order to discuss this accuracy improvement, we now com- 
pare method 3 to methods 1 and 2 using the SZ template used to 
create the mock data. The one dimensional parameter distributions 
obtained from the different methods are presented in figure [3] The 
black dotted lines represent the parameter distribution when we fit 
with pure primary CMB a signal that contains primary CMB alone. 
As expected, all 3 methods give unbiased parameters but their asso- 
ciated error bars (i.e. the accuracy of the parameter determination) 
differ. 

Unsurprisingly, adding an extra parameter (the SZ amplitude 
Asz) in the MCMC analysis with method 1 enlarges all the error 
bars (Fig. [3] long dashed blue curves) that would be obtained in 
a pure primary CMB analysis (black dotted line). Method 2 (dot- 
dashed green curves) with its explicit and strong dependency of the 
SZ amplitude on <r 8 improves the constraints on cr 8 , in comparison 
to the pure CMB analysis. At the same time, the accuracy on the 
spectral index n s is deteriorated. Increasing n s increases the CMB 
power at high t values but this can be slightly compensated by a 
decrease of the SZ amplitude through erg or Clbh. The solid black 
curves represent the expected precisions obtained with method 3. 
We note that they are improved for all parameters except n s as 
compared to the pure primary CMB analysis. The improvement is 
particularly important for <r 8 : a factor 2.5 better. The accuracy on 
Sli m h 2 and on the derived parameters f2 A , fl m , H is also improved 
by a factor 2. 

Method 3 provides better constraints on the cosmological pa- 
rameters because it allows to break degeneracies between models, 
as illustrated in figure [4^. This figure represents the 68.3% confi- 
dence region on parameters Q.i m h 2 and A s when the analysis is only 
performed on these 2 parameters. An analysis using only the pri- 
mary CMB signal gives the distribution represented by the red/dark 
ellipse. Using our SZ module to carry out an analysis on the SZ 



alone power spectrum (black ellipse) shows that A s and £la m h 2 are 
strongly degenerated but very differently with respect to a pure 
CMB, which means that the SZ contains some information com- 
plementary to that included in a pure CMB signal. The analysis of 
the CMB + SZ signals based on method 3 results in the two di- 
mensional distribution represented by the green/light ellipse that is 
smaller than the red one, taking advantage of the information in- 
cluded in the SZ power spectrum. 

Figure |4p represents two pure primary CMB spectra calcu- 
lated for the reference cosmological model (black solid line) and 
for a model in which the value of fidm^ 2 is the reference value plus 
one standard deviation that was obtained using method 3 (red dot- 
ted line). Since this standard deviation is much smaller (twice) than 
those obtained from a pure primary CMB analysis (see figure [3), 
these two CMB spectra are almost identical, showing the impossi- 
bility for a pure CMB analysis to distinguish between the two cases. 
The green continuous and dotted lines (figure[4p) respectively rep- 
resent the SZ power spectra for the reference and the shifted mod- 
els. The change in the £ia m h 2 value significantly affects the ampli- 
tude of the SZ spectrum. As a result, as shown by the blue curves, 
the CMB+SZ spectra are different and an analysis based on method 
3, that uses the information encoded in the SZ spectrum shape, al- 
lows to disentangle the two cosmological models and consequently 
better constrains the cosmological parameters than a pure primary 
CMB analysis. 

We note that the spectral index n s is the only parameter that 
is somewhat less precisely constrained through method 3 than it is 
in methods 1 or 2. An increase of its value raises both the CMB 
power and SZ spectra at small scales and lowers their power at 
large scales. However, while the CMB spectrum pivots at interme- 
diate scales, this happens at large multipoles for the SZ spectrum. 
As a result an increase of « s has two competing effects between 
these domains : it increases the CMB power and at the same time 
decreases the SZ power. The resulting CMB+SZ spectra for dif- 
ferent, but close, values of n s are thus more similar to each other 
than pure CMB spectra. That is why the spectral index is slightly 
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Figure 3. One dimensional parameter distribution. The coloured curves represent the distribution for the CMB+SZ signal fitted with the 3 different methods. 
Long dashed blue line : method 1, dot-dashed green line : method 2, solid black line : method 3. The SZ template used in method 1 and 2 is the same as that 
used to create the data. The black dotted line represents the parameter distributions when then signal contains only a pure primary CMB. The vertical red lines 
represent the input values of the parameters we used to create our mock data. 




Figure 4. Left panel: Two dimensional parameter distribution. The black ellipse represents the Qdm'' 2 / pairs obtained from the Markov chains analysis of 
the pure SZ spectrum. The orange (dark grey) ellipse represents the pairs obtained from the analysis of the pure CMB spectrum. The green (light grey) ellipse 
presents the pairs obtained using the CMB + SZ spectrum. Right panel: The black solid and green dashed lines are superimposed. They represent the pure 
CMB spectra respectively for the reference cosmological model and a one with a shifted value of Cl^ m h 2 (see text for details). The red thick solid and dashed 
lines represent the corresponding SZ spectra and the blue solid and dot-dashed line represent the CMB+SZ spectra. 



less accurately constrained in joint analysis of CMB+SZ than in an 
analysis of a pure CMB signal. Nevertheless this slight degradation 
of the accuracy on n s is negligible compared to the improvement of 
the precision of the other parameters. 

4.2 Effects of the intra-cluster gas physics 

As discussed in section[2], the computation of the TSZ power spec- 
trum involves assumptions about the gas physics. Either, on the one 
hand, hydrodynamical-numerical simulations are performed with 
given models for the gas evolution (adiabatic cooling, preheating, 
feedback, etc.), or on the other hand, theoretical computations of 
the power spectrum choose a model for the intra-cluster gas dis- 
tribution as well as scaling relations between the total mass and 
other cluster physical parameters (e.g. temperature). The large di- 



versity of the possible models describing the cluster gas properties 
reflects the difficulty to summarise in a simple parameterisation the 
complexity of the intra cluster gas physics. The cluster physical de- 
scription will inevitably affect the SZ spectrum and consequently 
the cosmological parameter estimation. 

In the following, we focus only on the theoretical computa- 
tions, for which modifying the cluster physics is straightforward, 
and we illustrate the effects of such modifications on the cosmolog- 
ical parameters. We restrict our study to one single intra-cluster gas 
distribution model: the y8-profile. We vary the value of the isother- 
mal gas temperature T e and the index /?, i.e. the overall amplitude 
and steepness of the density distribution respectively. We run an 
MCMC with method 3 to fit the data created with different values 
of p and T c . The resulting one-dimensional parameter distributions 
are presented in Fig. [5] 
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An electron distribution with /? = 0.8 (i.e. 20% larger than the 
reference value) corresponds to a more peaked individual SZ pro- 
file. This increases the overall SZ power and shifts the maximum 
of the SZ spectrum towards smaller scales. Using a primary CMB 
plus SZ spectrum calculated using the intra-cluster gas reference 
parameterisation to fit such data results in biasing the estimation 
of some cosmological parameters, namely Cl b h 2 (0.85 times the ex- 
pected accuracy (ea)), fldm/2 2 (2.6x ea) and n s (0.83x ea). This cor- 
responds to biases equal to 2.4 for C1 A and Q. m , 3.3 for cr 8 and 2.1 
for H in terms of the expected precision. A conservative 5% dif- 
ference on fi (fi = 0.7), still biases Adm'? 2 value at the level of 0.5 
times the expected accuracy. The values of S1 A and Q m , erg and H 
are also biased at several tenth of the expected precision (0.4, 0.8 
and 0.4). 

Overestimating the electronic temperature at a 7% level (accu- 
racy on the electron temperature obtained from X-SZ clusters, see 
iReese et al., 2002; Bo namente et alJuOOfj) in the data analysis re- 
sults also in a biased estimation of some cosmological parameters, 
namely Q.\,h 2 (0.6x ea), CIdmIi 2 (1.7x ea) and n s (0.5x ea). This 
translates to biases of 1.5 for f2 A and Si m , 2.2 for cr 8 and 1.4 for Hq. 



5 CONCLUSIONS 

In this paper we investigate different methods to jointly fit the pri- 
mary CMB and SZ signals. We emphasize that methods using a 
power spectrum template to describe the SZ contribution are likely 
to bias the cosmological parameter values in the case of an inap- 
propriate choice of the SZ template. This is due to the frozen shape 
of the SZ spectrum that is inconsistent with changes in the cosmol- 
ogy. We also show that a joint CMB-SZ analysis with a full cosmo- 
logical dependency of the SZ spectrum does not suffer from such 
biases. Moreover it improves the estimate of cr 8 and Cli m h 2 (2.5 and 
2 times respectively) with respect to a pure CMB analysis. In that 
respect, such a method can be presented as a coherent analysis of 
the primary CMB plus SZ signal. However, we point out that our 
incomplete understanding of the intra-cluster gas distribution and 
properties can result in errors in the calculation of the SZ angular 
power spectrum and is likely to bias the cosmological parameter 
estimation up to 2 sigmas for fl m and H and even more for <x 8 . 
This reinforces the need for better constraints on the description of 
the galaxy cluster gas properties. 
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Figure 5. The curves represent the one dimensional parameter distribution for the CMB+SZ signal fitted with method 3. The solid black line represents the 
parameter distribution when the gas description used in the fitting method is in agreement with the one employed to create the data. The green long-dashed 
(dot-dashed) line presents the case in which the gas description used to create the data uses a /? = 0.8 (0.7) parameterisation. The blue dotted line stands for 
the case where the electronic temperature is overestimated at a 7% level in the intra cluster gas description used to fit the data. The vertical red lines represent 
the input values of the parameters we used to create our mock power spectra. 
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